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ADAPTIVE DETECTOR ARRAYS FOR OPTICAL COMMUNICATIONS RECEIVERS 

CROSS-REFERENCE TO RELATED APPLICATION(S) 

10 This application claims the benefit of U.S. provisional application number 60/253,61 0 filed 

on November 28, 2000, the content of which is incorporated herein by reference. 

FIELD OF THE INVENTION 

This invention is directed to improved optical receivers, and particularly to high data-rate 
15 optical receivers utilizing wide-band optical detector arrays capable of correcting for signal 
degradation caused by atmospheric turbulence. 

O 

f BACKGROUND OF THE INVENTION 

P It is well known that ground-based reception of optical signals suffers from degradation of 

20 hi the optical phase-front caused by atmospheric turbulence. This turbulence leads to a reduction in the 
Ut e ff ec tive diameter of the receiving telescope, and to random fluctuations of the receiver's "point 
f=A spread function" (PSF) in the focal plane. 

hi For example, the diffraction-limited field of view (FOV) of a receiving telescope can be 

I*| taken to be approximately 0 d i = A/D R , which, for a 3-m aperture and 1 um wavelength, translates to 
25 ¥ s 0.33 farad. If the effective focal length of the telescope is 6 m (implying an F/2 instrument), then a 
diffraction-limited PSF of 2 urn diameter, or 0.002 mm, will be produced in the focal plane. Thus, 
under ideal conditions a very small detector could be used to collect virtually all of the signal energy, 
while at the same time spatially filtering out most of the background radiation. 

However, atmospheric conditions rarely permit diffraction-limited operation of large 
30 telescopes; even under "good" nighttime seeing conditions, the phase of the received signal field 
tends to become uncorrelated over distances greater than 20 cm, deteriorating to as little as 2 to 4 cm 
during the day. Under these conditions, the dimensions of the PSF in the focal-plane tends to 
increase inversely with coherence length, as if the diffraction-limited telescope were correspondingly 
reduced; the telescope still collects all of the signal energy propagating through its physical aperture, 
35 but the collected signal energy is redistributed into a much larger spot in the focal plane. In 
conventional receivers, the receiver's FOV is increased proportionally to collect the signal. 
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5 However, this increase in the receiver's FOV leads to a corresponding increase in the amount of 

interfering background radiation offsetting much of the performance gain. 

Some attempts have been made to utilize signal-processing hardware to reduce the 

deleterious effects of atmospheric turbulence on receiver performance. However, due to the 

limitations of the chosen algorithms and the electronics utilized, only the average PSF of the 
10 received signals over a relatively long time period have been successfully processed. Although 

some performance improvement has been seen from receivers utilizing these time-average PSF 

signal-processing techniques, because the PSF can change on the order of milliseconds much of the 

detailed information from the processed transmissions is lost. 

Accordingly, there is a need for an optical communications receiver capable of correcting the 
15 signal degradation from atmospheric turbulence instantaneously and without significantly increasing 

interference from background radiation. 

•F3fs? 

:f% 

% SUMMARY OF THE INVENTION 

s w 

1*1 The present invention is directed to an optical communications receiver comprising a wide- 

•i % % 

20 s band optical detector array and a high-speed digital signal processor assembly programmed to 
jf! operate on the raw data from the detector array in real-time to ameliorate the effects of atmospheric 
111 turbulence on the performance of the optical receiver while operating within the terrestrial 
'm atmosphere, or while attempting to communicate through any similar turbulent medium. 
?a * In one embodiment, the detector array is designed for communications applications such that 

25 the signal processing assembly processes an array of parallel outputs, one from each detector 
element, instead of the usual serial row or column readouts typical in imaging applications. In such 
an embodiment, each detector element preferably has a bandwidth commensurate with the data 
requirements. Any parallel output detector array may be utilized in such an embodiment, such as, for 
example, an array of photo-multiplier tubes (PMTs) capable of counting individual photons or an 
30 array of "avalanche photodiode" (APD) detectors with individually accessible outputs. 

In one embodiment, the signal processing assembly samples each output from the detector 
array at the Nyquist rate, and uses these samples to measure the average signal and background noise 
energies over each detector element. In such an embodiment, since the coherence time for the 
turbulence-degraded fields is on the order of several milliseconds, the integration time for this 
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measurement can last up to a millisecond or more, during which time a great many data-symbols are 
received. 

In the above embodiments, several possible algorithms can be applied via the signal 
processing assembly to combine the detector outputs advantageously, i.e., to maximize the 
probability of correct detection. For example, in one embodiment, a logarithmic function of the ratio 
of signal and noise energies is applied to each output for a short time corresponding to the 
coherence-time of the channel. In an alternative embodiment, a time-varying "mask" is computed 
and applied to the array such that only those elements of the array with the greatest signal content are 
utilized. 

In still yet another embodiment, the invention is directed to a method of ameliorating the 
effects of atmospheric turbulence on the performance an optical receiver utilizing the detector array 
and signal processing assembly described herein. 

BRIEF DESCRIPTION OF DRAWINGS 

These and other features and advantages of the present invention will be better understood by 
reference to the following detailed description when considered in conjunction with the 
accompanying drawings wherein: 

Figure la is a representation of a receiver's point spread function under ideal conditions; 

Figure lb is a representation of a receiver's point spread function under atmospheric 
turbulence conditions; 

Figure 2 is a schematic of an optical receiver according to the present invention; 

Figure 3 is a flowchart of an embodiment of a signal processing method utilized in the optical 
receiver according to the present invention; 

Figure 4 is a representation of a receiver's point spread function overlaid with a signal mask 
by a signal-processing assembly according to the present invention; 

Figure 5 is a graphical representation of the performance of an optical receiver according to 
the present invention in terms of average number of received photons; 

Figure 6 is a graphical representation of the performance of an optical receiver according to 
the present invention in terms of average number of received photons; 

Figure 7a is a graphical representation of the performance of an optical receiver according to 
the present invention in terms of photon efficiency; 
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5 Figure 7b is a graphical representation of the performance of an optical receiver according to 

the present invention in terms of photon efficiency; 

Figure 8 is a graphical representation of the performance of an optical receiver according to 
the present invention in terms of number of ranked detector elements included; 

Figure 9a is a graphical representation of the performance of an optical receiver according to 
10 the present invention with various PPM dimensions, and background noise levels; 

Figure 9b is a graphical representation of the performance of an optical receiver according to 
the present invention with various PPM dimensions, and background noise levels; 

Figure 9c is a graphical representation of the performance of an optical receiver according to 
the present invention with various PPM dimensions, and background noise levels; 
15 Figure 9d is a graphical representation of the performance of an optical receiver according to 

! % the present invention with various PPM dimensions, and background noise levels; 

ttl Figure 10 is a graphical representation of the performance of an optical receiver according to 

?y 

HI the present invention for different realizations of the instantaneous point spread functions; 
\ Figure 1 1 is a graphical representation of the performance of three different signal-to-noise 

20 Iti ratio measures considered in an optical receiver according to the present invention with; 

!„! Figure 12 is a graphical representation of the performance using different signal-to-noise 

J! ratio measures according to the present invention; and 

III Figure 13 is a graphical representation of the effect of real-time signal-to-noise ratio 

13 

y. estimations on the performance of an optical receiver according to the present invention. 

25 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention is directed to an optical communications receiver comprising a wide- 
band optical detector array and a high-speed digital signal processor assembly programmed to 
operate on the raw data from the detector array to ameliorate the effects of atmospheric turbulence 

30 on the performance of the optical receiver in real-time while operating within the terrestrial 
atmosphere, or while attempting to communicate through any similar turbulent medium. 

An example of the increase in the effective dimensions of an exemplary receiver point-spread 
function over its diffraction-limited value as a result of atmospheric turbulence is shown by 
comparison in Figures la and lb. The hypothetical signal shown in these figures correspond to a 

35 telescope having aim aperture, an 0.3 m central obstruction, and a 4 cm atmospheric coherence. 
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5 As shown in Figure la, under ideal conditions the undistorted signal generates a diffraction limited 
point spread function (PSF) in the detector array. However, as shown in Figure lb, under the 
influences of atmospheric turbulence, the point-spread function is degraded and enlarged. 

Under such conditions, in order to collect all of the signal energy, the dimensions of a single 
optical detector must be made large enough to encompass the degraded point-spread function as well 
10 as its random excursions in the focal plane, which tend to change on time scale of 10 to 100 ms. 
Thus, the active area of the detector must be made large enough to encompass most of the signal 
energy most of the time. However, a large detector implies a large receiver field of view, which in 
turn implies a corresponding increase in the amount of background radiation admitted into the 
receiver. That, in turn, degrades communication performance. 
15 m Although in the above example a telescope of 1 m diameter with a secondary obstruction of 

*lS 0.3 m has been assumed, it should be understood that similar results hold for larger telescope 
yy diameters as well as long as the focal-plane signal distribution is dominated by turbulence effects. 
Jil Since the diffraction-limited PSF is inversely proportional to the telescope diameter, this condition 
W generally will be satisfied for larger telescope diameters; however, the background and signal 
20 s energies must be properly scaled to account for the larger collecting area. If the receiver optics are 

i..H. 

\2 not diffraction limited but have significant imperfections, generating a PSF with dimension that are 
!'W comparable to the turbulence-induce distribution, then more detailed modeling of the signal 
Ft distribution becomes necessary. 

■ " A conceptual block diagram of an optical photon-counting array receiver according to the 

25 present invention is shown in Figure 2. The receiver 10 consists of a collecting aperture 12 and 
optics 14 to focus the collected fields onto the focal plane 16 of a detector array 18 comprising a 
plurality of detector elements 20, such as individual photomultiplier tubes or avalanche photodiode 
detectors, which respond to the impinging fields. The individual output voltage 22 from every 
detector element 20 of the array 18 are converted to numeric samples, which then are operated on by 
30 a signal-processing assembly 24 that performs the required signal processing algorithms in order to 
optimize the receiver performance in real-time. 

In the above embodiments, any collecting aperture 12 capable of collecting optical signals 
from an external source and transmitting those collected signals to the receiver optics 12 may be 
utilized in the current invention. For example, 10 cm telescopes for terrestrial repeaters, 5 to 10 m 
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5 communication apertures, or 1 m aperture optical communication telescope laboratory (OCTL) 
receivers. 

Likewise, any optics package 14 suitable for focusing the collected fields onto the focal plane 
16 of the detector array 18 may be utilized in the current invention. 

Although the above embodiment only describes the use of photomultiplier and avalanche 
10 photodiode detector elements 20, any suitable optical detector may be utilized in the current 
invention. For example, even PIN diode arrays could be utilized when adequate signal energy is 
available. 

The signal processing assembly 24 of the current invention may comprise any electronic 
components suitable to allow signal analysis to occur at sufficiently high-speed to allow 

15 instantaneous or real-time analysis of the PSF, i.e., which allows signal processing of each detector 

Q 

yT| element output at the Nyquist rate not at a signal averaged rate. For example, high-speed digital 
"*J signal processors (DSPs), FPGAs, or dedicated ASIC designs could be utilized. 
IP The signal-processing assembly 24 optimizes the receiver performance by selecting which 

l s i elements 20 of the detector array 18 should be used to collect data at any arbitrary instantaneous 
20 -"^ point of time. Data from detector elements 20 with signal content sufficiently over the background 
\dk signal level to improve the signal statistics are then combined and data from elements 20 with signal 

m levels at or near the background signal level are rejected and the result is output by the signal- 

i if 

'ff processor assembly 24 to a decoder 26 which detects the incoming signal at that arbitrary time point 
from the signal-processor assembly output. 

25 Figure 3 depicts a flow-chart of an embodiment of the process utilized by an exemplary 

signal-processor assembly 24 according to the current invention to optimize the receiver 
performance. During operation, the signal-processor assembly 24 first obtains a vector array of data 
outputs from the detector array 18 at Step 1 and then at Step 2 operates to sample the data to measure 
the background and signal noise energies and then compute an estimate of the signal intensity and a 

30 data weighting parameter for each output based on one of several possible processing algorithms to 
increase the probability of correct signal detection. At Step 3, the detector outputs are multiplied by 
their corresponding calculated data weighting parameter and outputs the weighted data. At Step 4 
the combined weighted data is utilized to detect the transmitted signal and the signal is output to the 
user at Step 5. 
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5 Any data weighting algorithm and parameter may be utilized to optimize the signal received 

by the processor 24 from the detector array 18, such as, for example, applying a logarithmic function 
of the ratio of signal and noise energies to each output for a short time corresponding to the 
coherence-time of the channel, or alternatively, computing and applying a time-varying "mask" to 
the array that effectively uses only those elements of the array with the greatest signal content. 

10 These exemplary algorithms and the signal-processing method utilized in the current invention, and 
schematically shown in Figure 3 are described in greater detail below. 

In a conventional single-detector element receiver, the receiver measures the number of 
photons contained in a received field by producing a stream of free electrons at its output terminal in 
response to the absorbed photons. If the occurrence time of each pulse can be measured, and if the 

15 q amplitude of each pulse is normalized to unity, then the count record consists of positive integer- 
^ valued jumps occurring each time a photon is detected. Therefore, the count record or count 

sJ3 accumulator function is a monotonically increasing function that contains all of the information 

01 

|i| present in the detection process and the conditional density of the total number of counts, N, given 
7^ an intensity function X(t), can be expressed as: 

20 *■ 

fly p[N(t)\A(t);0 <t<T] = exp[- )x(t)dt\ (la) 

13 

where N(T) = 0, and 

25 p[N(t)\A(t);0<t <T] = 

where N(T) > 1. Where the set {wj represents the occurrence times of the detected pulses; N(t) is 
called the count accumulator function of the process over the time interval (0,T); and X(t) = 
jJ" A dxdyA.(x,y;t), where A is the detector area and X(x,y;t) is the intensity process over space and 
30 time. 

However, this expression only describes the output of a single detector element responding to 
an optical intensity. When an array of detectors is used to detect the optical field, Equations la and 



i=0 



exp 



jk(t)dt 



(lb) 
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5 lb must be generalized to enable an unambiguous description of the output of each detector 
element 20. 

Consider a rectangular array of detector elements 18 of Kx L detector elements 20, as shown 
in Figure 2. For some applications, such as finding the center of the signal intensity distribution, it is 
important to know the location of each detector element within the array. In such an embodiment, 
10 the subscripts mn, where 1 < m < K and 1 < n < L to denote the position of the detector element 
within the array. Thus, the sample function density defined in Equation 1 can be written as: 

pKMKM0^t<T] (2) 



15 which now represents the output of a particular element of the array. Note in such an embodiment, 
=35 ^mn(t) can be viewed as that portion of a spatially distributed intensity function intercepted by the 
^ mnth detector element according to: 

U Kn(0= \\dxdyX{x,y;t) (3) 

si- 

20 y, 

where A mn is its effective area. Note that if the spatial intensity distribution is known, and the 
||l location and size of each detector element also are known, then conditioning on the spatial intensity 
distribution is equivalent to conditioning on the array of intensity components, each of which is still 
a function of time. Assuming that each array element observes the sum of a signal field plus 
25 multimode Gaussian noise field with average noise count per mode much less than one, the array 
outputs can be modeled as conditionally independent Poisson processes, conditioned on the average 
signal intensity over each detector element. Hence, the joint conditional sample function density of 
the array can be denoted as: 

30 p[N(t)\Z{t)$ <t<T] = f[fl p[N m „ (t)\X mn (0;0 < t < T] (4) 

m=\ n=l 



where N(t) = (Nn(t),Ni2(t), . . ., NklCO) and each component on the right-hand side of Equation 4 is 
defined in Equation 1 . 
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5 Accordingly, a signal propagating through the atmosphere and through the receiving 

optics 12 is transformed into a space-time intensity function in the detector plane. The receiver 10 
also collects background energy from all directions, which is assumed to contribute an additional 
constant intensity, Xb per detector element. The integrated intensity is then given by Equation 3, 
above plus the additional constant background level. 

10 Once the total intensity is measured by the detector the signal-processor operates on the 

signal to instantaneously determine the information contained therein. To accomplish this signal 
processing, in one embodiment, at the end of T sec, the signal processor assembly estimates the 
signal intensity and then computes the probability of having received the observed array of count 
accumulator functions and selects that message corresponding to the greatest probability of having 

15 been received. The spatial component of the process typically consists of calculating the greatest 
O log-likelihood function, A(T), having a spatial portion given by an equation of the form: 



firs K L 



s.mn 



m=\N=\ {(observational interval) 



1 + - 



N mn (5) 



lis 



20 |4 where is the total number of photons occurring over the mnth detector element. It should be 
pjJi noted that the log-likelihood function of Equation 5, above, mathematically relates only the spatial 
IS portion of the function and that such a function may require an additional modulation dependent 
=u term. In addition, the spatial intensity must be estimated often enough to capture any fluctuations in 
the modulation, but must be averaged over sufficient time to ensure adequate signal statistics (in 
25 most cases -msec ). Rewriting the logarithmic functions or weights in Equation 5 as u mn , allows the 
log-likelihood function to be rewritten as: 



K L 



A(T) = ^u mn N ma (6) 



30 In this form, it is clear that the log-likelihood function is composed of sums of a random 

number of weights from each detector element, for example, the mrath detector element contributes 
an integer number of its own weight to the sum. 
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5 According to the preceding analysis, in a receiver system according to this invention, 

elements of the detector array containing much more background than signal intensity do not 
contribute significantly to the error probability, since the output of these detector elements is 
multiplied by weights that are close to zero. Accordingly, in an alternative embodiment of the signal 
processing assembly 24 the processor first ranks the detector elements, starting with the one 

10 containing the most signal energy and followed by every other detector ordered according to 
decreasing signal intensity. The processor then computes the probability of error for the first 
detector element plus background; then for the sum of the signal energies from the first two detector 
elements (plus background for the two detector elements), and so on, until the minimum error 
probability is reached. In such an embodiment, each set of detectors may be considered to be a 

15 single detector, so that no weighting is applied to account for variations in the signal distribution 

CJ over the detector elements included in that set. The set of detector elements that achieves the 

i.fi 

% minimum probability of error is the best synthesized single detector matched to the signal-intensity 
■*| distribution and is the one chosen as the signal. 

U-'r 

III However, it should be understood that the above system requires intensive calculations. 

20 sT? Accordingly, in one alternative embodiment, the logarithmic weights are partitioned into two 
f classes: large weights are assigned a value of "1" and small weights are assigned a value "0." 
I** Effectively, such a signal-processing method creates a "signal mask" that groups a set of the highest 
•jj^ ranked detector elements together into a synthesized single detector as shown schematically in 
P Figure 4. In turn this synthesized detector changes with the changing channel. 
25 Although only a binary approximation has been discussed above, any other suitable 

approximation could be applied to the logarithmic rates. For example, the computation of log(l+x) 
could be replaced by the computation of "x" in some cases, namely when the signal-to-noise ratio 
over each detector element is small. 

Because computing the error probability for each increasing subarray requires a great deal of 
30 computation time, particularly when large detector arrays are used, in an alternative embodiment of 
the invention, the calculation takes into account the number of detectors and the total average signal 
and background energies. Utilizing such a method three different computations analogous to signal- 
to-noise ratios (SNR) were created according to: 
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SNR 2 (/) = fei^pj (8) 



and 



T 



™ 3 (/) = J ^ '-" (9) 

10 

^ where the index / represents the number of detector elements over which these functions are 
tff maximized, 1 < / < 256 for the current example. 

sQ Although the above signal processing method and apparatus may be utilized with any 

Si intensity modulation technique, such as for example, binary PPM and on-off keying, in an 
15 p| exemplary embodiment the technique will be modeled on a M-ary PPM technique. 

% In M-ary PPM modulation, a signal pulse of duration x sec is transmitted in one of M 

y, consecutive time slots, resulting in a PPM symbol of duration T = xM sec. Under the hypothesis that 
-Ljf the signal pulse is contained in the z'th time slot, the signal intensity is given by: 



20 *«(/) = „ \ (10) 

0 : else 



As before, after propagating through the atmosphere and through the receiving optics, this 
temporal intensity function is transformed into a space-time intensity function in the detector plane. 
The integrated intensity over the mnth detector according to Equation 3 can then be designated as 
25 hf !) mn(t) in order to incorporate the hypothesis dependence. Where the integrated intensity can be 
written as: 

IX „» (0 + K '■ (i -l)r<t< it 
^o (0 = J * > (11) 

I k b : else 
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5 In this exemplary embodiment, it is assumed that each of the M messages is equally likely to 

be transmitted with probability M 1 , and that each message generates a unique vector of detector 
array intensities at the receiver, denoted by X (i) (t) = (X & n (t), X Q \ 2 (t), . . . , >- (i) KL(t)). As before, at the 
end of T sec, the signal-processing assembly computes the probability of having received the 
observed array of count accumulator functions and selects the message corresponding to the greatest 

10 probability of having been received. Equivalently, the decoder selects the message corresponding to 
the greatest log-likelihood function, K(T), evaluated after T sec and conditioned upon the signal 
occurring at the ith time-slot: 



i0 



a,, (r) = \n{ P [N(tp\t);0 <t < rj = EE ln {^k«(0|^i(0;0 < t < t} 

m=\ n=\ 

(12) 



K L ( ir 



= JZ - J ^(0&+ Y^^niWj,™) + {terms that depend only onX b ) 



m=l n=l 

TU 
15 W 

III 

s where Wj_ mn is the occurrence time of the 7'th photon over the mnth. detector element within the same 
fT time slot. 

111 Assuming constant signal intensity over the ith slot, independent of the value /, and ignoring 

fi the terms that depend only on Xb (which hence do not convey any information about the 
20 ^ transmission), the log-likelihood function reduces to: 



A,(70 = yElnl + ^V;£ (13) 



m=l n=\ 



V X b J 



where vP mn is defined as the total number of photons occurring over the mnth detector element 
25 during the ith time slot. Note that in this exemplary embodiment, with constant signal intensities the 
actual arrival times of photons within each time slot do not contribute to the decision; hence, only 
the total number of detected photons, f/^ mn , matters. Given that the intensity over each detector 
element is known, the ith log-likelihood function consists of the sum of a logarithmic function of the 
ratio of signal and background intensities from all detector elements over the ith pulse interval, 
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5 multiplied by the total number of detected photons; the optimum detection strategy is then to select 
the symbol corresponding to the greatest log-likelihood function. 

In such an M-ary PPM system, two distinct cases arise, when signal is present and when 
signal is not present. When signal is present then the intensity over the mnth. element of the detector 
array is equal to the sum of the signal intensity X mn and the background intensity Xb, whereas when 
10 signal is not present the intensity over each detector element is simply equal to the background kb. 

In such a M-ary PPM system, the received information is processed and decoded correctly if 
the weighted sum for the signal slot (q) exceeds the weighted sum for every other , i.e., non-signal 
slot, such that: 



15 

y 



Am>A ; (r) 



(14) 



20 s 



25 



However, since the weights for each log-likelihood function are the same, it is possible that the 
maximum sum of weights occurs over two or more times, one of which is the true signal time, 
creating a "tie." In the case of a "tie" for the largest weighted sum in such a system, random choice 
is utilized among the largest likelihood functions: if there are r nonsignal slots tied for biggest with 
the correct signal slot, then the probability of selecting the correct symbol in the presence of (r + 1) 
ties for biggest is {r + l)" 1 . However, there are, (M - 1) taken r at a time, ways these ties can occur 
among M slots. On the other hand, if none of the detectors registers a photon over any of the M 
slots, then we have to make a random choice among M possibilities, which yields a correct decision 
with probability M 1 . Assuming equiprobable signals, the probability of correctly decoding the 
received signal is denoted in Equation 15, below. 



30 



P M {C) = P M {C\H q ) = 



Af-l 

I 



f 1 "\ 



r + lj 



M-l 



' k=l 



l*q 



M-l-r 



> + M^\p q (a 0 \H q \ Pi (a 0 \Hj- 1 ' : 



(15) 



where p q and p t , i * q, refer to the probability densities corresponding to the signal and non-signal 
hypotheses, respectively. 



13 



C766/JWP/47456 



5 Alternatively, counting all "ties" as errors yields a tight lower bound on the probability of 

correct detection, and yields an easier to calculate probability function: 



P M (C)>PUQ = tpM 



H 



4=1 



+ 



^IpMHjp^Hj' 1 } (16) 



10 In turn this lower bound on the probability of correct detection yield the following upper bound on 
the error probability, which may be utilized to evaluate a variety of proposed detector array schemes: 



P»{E) = \-Pl{C)>P M {E) 



(17) 



15 J3 

m 
m 

l.i 



20 ill 



Although the above analysis has assumed an optimal detector regime, such a M-ary PPM 
system might also use a suboptimal or adaptive synthesized single detector, as described above. In 
such an embodiment, the probability of correct detection can be obtained directly from Equation 15, 
by setting a* = k in the probability densities and assuming constant signal and background intensities 
over each time slot, yielding: 



k\ 



-(V+4t ) 



and 



(18) 



(19) 



25 Direct substitution of these Poisson densities into Equation 15 yields a probability of correct 

detection according to: 



M-l 

y 



' 4=1 



4-1 



M-l-r 



(20) 
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5 

However, when the array contains a large number of detector elements, the computation of 
the probability density of the weighted sum of Poisson random variables by the signal processing 
assembly may become prohibitive. In such a case, an alternative approximating signal processing 
method may be utilized to substitute the true discrete density with a simpler continuous density. For 
10 example, in one alternative exemplary embodiment, a Gaussian approximation to the discrete density 
of the weighted sum of Poisson random variables is derived from the characteristic function of the 
discrete density. 

The various signal processing techniques detailed above were modeled utilizing an 
exemplary M-ary pulse-position modulation (PPM) data transmission protocol. In this protocol one 
15 of M intensity functions is received by the receiver 10 and the receiver attempts to determine the 
,JJ correct symbol based on observations of the array of count accumulator functions over each of M 
*f time slots. It is assumed for the purpose of this exemplary embodiment that the symbol boundaries 
§1 are known and that the arrival time of each detected photon and total number of detected photons 

z 

i*j can be stored for a limited duration of time necessary for processing by the signal-processing 
20 W assembly 24. 

14, First, a comparison of the performance of signal-processing assemblies utilizing the 

optimally weighted array receiver and the adaptive synthesized single-detector signal-processing 
Ill methods described above, were carried out for two model receivers exposed to average background 
y, energies of K b = Xtfi = 0.1 and 1.0. 

25 Two different signal models were used: 1) a simple test model wherein only 5 of the 16x16 = 

256 total detector elements were assumed to contain signal energy while the rest were assumed to 
contain no signal; and 2) a more realistic 16x16 detector array model wherein the signal distribution 
over the array was simulated using a Kolmogorov turbulence model in which all 256 detector 
elements potentially contained some signal. 

30 For the test model, the proportions of the total average absorbed signal energy, K s = X s x, over 

the five detector elements were assumed to be (1.0, 0.3, 0.2, 0.05, 0.02). The optimally weighted 
array receiver (Equation 16) and the synthesized single-detector receiver (Equation 20) signal 
processing methods were then evaluated with this model and compared with results obtained via 
Monte Carlo simulations. The results are shown in Figure 5, as a function of the total average 

35 absorbed signal energy, K s . From this plot, it is evident that optimal weighting yields somewhat 
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better performance than the suboptimum single-synthesized or signal mask array, and that greater 
improvements occur at greater background intensities. However, the improvements due to the more 
complicated signal-processing utilized in the optimally weighted array are only about 0.3 dB at an 
error probability of 0.001 for the high background case. 

In Figure 6, a realistic spatial distribution of the signal intensity over the focal plane was 
generated using Kolmogorov phase screens. Monte Carlo simulations were performed to evaluate 
the error probability for the optimally weighted array. The Gaussian approximation to the error 
probability defined by: 



K L 



71 = > > U X X 



(21) 
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and 
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was also evaluated for the binary PPM case, M= 2, with average background energies (per detector 
element) of K b = % b x = 0.1, 1.0 and 5.0 photons/time slot, as a function of the total average signal 
energy. 

From this plot, it can be seen that the Gaussian approximation is close to the exact values 
obtained from Monte Carlo simulation and that good agreement is obtained even for small 
background energies per detector element, as direct comparison with the simulation results indicates. 
In fact, the Gaussian approximation embodiment provides useful results over the entire range of 
background and signal energies represented in Figure 6. The performance of the synthesized signal 
masked process is also shown; it can be seen that its performance is well within tolerable limits. 

Both analytical calculations and Monte Carlo simulations were performed in order to obtain 
PPM error probabilities for the receiver according to the invention. Performance of the optimally 
weighted array receiver embodiment was obtained from simulations; for each PPM symbol, M 
Poisson random variables with the proper statistics were generated, the optimum weights were 
applied, and the symbol corresponding to the largest observable selected. Simulated turbulence- 
degraded signal distributions were then generated over the 16x16 detector array for all subsequent 
results. With no loss in generality, the transmitted symbol was always assumed to be the one 
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5 corresponding to a signal pulse in the first slot. The detector process was repeated a large number of 
times (until 100 errors were accumulated) and repeated for increasing average signal energy with 
various background levels. The results of this test are shown in Figure 7a and 7b for M = 2, 16 and 
256. The probability of bit error is shown as a function of the receiver's photon efficiency, p, which 
is a measure of the average number of bits of information carried by each photon. It can be seen 
10 from the plot that, with background levels from 0.1 to 1.0 photons per slot, p from 0.3 to 0.5 
bits/photon can be achieved with 256 PPM signaling at uncoded symbol-error probabilities of 
around 0.001 to 0.01. 

In order to generate a spatial distribution of the signal incident upon the detector plane, a 
sample field was generated using a Kolmogorov phase-screen program, resulting in a matrix of 
15 complex signal amplitudes. For the simulations, an atmospheric correlation length of r Q = 4 cm was 
W assumed, which implies that the results should apply to any receiving aperture that is much greater 
*P than this correlation length. The field intensity generated in the detector plane by the simulation 
J then was integrated over the elements of a 16x16 detector array, which was assumed to encompass 
[W the extent of the signal distribution in the detector plane. The detector signal intensities are 
20 Id normalized so that for the mnth detector an average number of absorbed signal photon of X s _ mn x is 
m obtained. A constant average background photon energy of X b i is assumed over each detector 

^ element. 

i .if 

g| For a given sample function of the intensity distribution, the 256 detector elements were 

P sorted in decreasing order of average signal energy, and Mary PPM symbol-error probabilities were 

25 calculated for increasing numbers of detectors, starting with the first detector, using the synthesized 
single-detector signal processing embodiment of the current invention. Figure 8 shows the symbol- 
error probability for binary PPM, M = 2, as a function of the number of detector elements used for 
the case K s = 10 and K b = 0.1 (that is, the total average signal photons absorbed by the entire array is 
10 and the average number of background photons per detector element is 0.1.) It can be seen from 

30 the plot that, for this case, the smallest error probability of 0.0049 is achieved by assigning unity 
weight to the first 1 1 detector elements containing the greatest signal intensities and zero to all the 
rest. 

In Figures 9a and 9b, binary PPM symbol-error probabilities are shown as a function of total 
average number of absorbed signal photons for four cases: 1) when the optimum number of 
35 unweighted detector elements is used, 2) when the optimally weighted array is simulated, 3) when 
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5 all 256 detector elements are given unity weight (synthesizing a large, nonadaptive single-detector 
element), and 4) when an ideal adaptive optics system succeeds in concentrating all of the available 
signal energy into a single detector element, which then is the only detector element that is observed. 

Using the same focal-plane signal distribution as before, error probabilities were computed 
for average background photon counts of 0.1 and 1.0, shown in Figures 9a and 9b, indicating 
10 performance gains by the adaptive detector array receiver according to the current invention over a 
single large nonadaptive detector of 2 and 2.8, respectively, at an error probability of 0.001, 
corresponding to 3 and 4.5 dB of performance improvement. Meanwhile, when compared with the 
ideal adaptive optics receiver that concentrates all of the collected signal energy in a single element 
of the array, the gains are 3.8 and 8.2, corresponding to 5.9 and 9.1 dB of improvement. However, 
15 between tne optimally weighted array embodiment and the synthesized single detector or signal 
■J% mask embodiment of the current invention, there is only a 0.3 dB improvement at a symbol error 

probability of 0.001, even with a relatively high background energy of K b = 1.0. 
0] Similar gains are evident in Figures 9c and 9d, which represent the symbol-error probability, 

hj P(SE), of the optimized subarray observing 16-dimensional (M = 16) PPM. The accuracy of the 
20 ^ union bound evaluated for the case X b x » 1 is evident, especially at the lower error probabilities, 
f * As shown in Figure 10, performance improvements are also observed for several different 

HJ focal-plane distributions at an average background energy of one photon per detector per slot to 
S verify that the above results were typical. This plot shows that three out of four simulations yielded 
t«* performance comparable to that shown in Figure 9b, requiring approximately 26 signal photons to 
25 achieve an error probability of 0.001 while utilizing 9 to 13 elements of the array in the region of 
error probabilities examined. Accordingly, receiver performance is independent of array distribution 
under the same atmospheric and receiver parameters. 

Figure 11 shows the three functions versus / for the case X s x = 10 and X b x = 0.1, assuming 
Poisson statistics. As can be seen from the plot, the maximum values for these three measures are 
30 reached when / equals 10, 22, and 7 detector elements, resulting in binary PPM symbol-error 
probabilities or 5.978 x 10" 3 , 5.84 x 10" 3 , and 1.049 x 10 " 2 , respectively, as compared with the 
performance of the true optimum subarray of 4.884 x 10 " 3 calculated from the actual error 
probabilities and achieved with 15 detector elements. 

Figure 12 in turn shows how subarrays obtained using the three different SNR measures 
35 given by Equations 7, 8 and 9, respectively, compare with the performance of the true optimized 
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5 subarray over a range of signal energies for a particular signal distribution. As can be observed from 
the plot, both SNRj and SNR 2 yield performances comparable to that of the optimized subarray and, 
therefore, could be used to optimize the subarray dimensions in real time. 

Finally, thus far the above evaluations of the performance of the receivers according to the 
present invention were obtained under the assumption that the true value of the average signal and 
10 background photons absorbed by each detector element are known. Accordingly, in these tests the 
sorting of the detector elements is based on the true signal energies. However, in many systems the 
signal energies will change with time due to turbulence. Therefore, the case when the signal 
energies were not known a priori, but had to be estimated was tested. The results of the simulations 
in which actual Poisson deviates were generated for each array element, and the mean signal 
15 energies estimated from the observed outputs, are presented in Figure 13 for binary PPM. 

U For each detector array element, Poisson random variable were generated for the Mary signal 

CI and background slots with average intensities obtained from the Kolmogorov phase-screen output 
|| plus a specified level of background light. These deviates were generated for L PPM symbols and 

flJ the slot outputs added for each detector element, resulting inmxn statistics according to: 

m 

20 |_ E j 

S L M-l 

f* L=ZI^(',;) (23) 

■ 

where X mn (ij) is the output of they'th slot of the z'th PPM symbol. These statistics were sorted as 
before from largest to smallest. The average number of signal photons was estimated from these 
25 statistics as: 



A s>m „r = max 



(24) 



where it is assumed in this exemplary embodiment that the actual background intensity can be 
estimated accurately, both because it is essentially constant and because typically there is significant 
dead time between PPM symbols to allow for transmitter laser recovery, which can be used to 
estimate the background intensity accurately since no signal photon are present. During operation 
once the average signal and background energies have been estimated, SNRfl) may again be 
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5 maximized over the number of elements, /, as before in order to obtain the optimum detector 
subarray. 

Results for binary PPM with estimated signal energies are shown in Figure 13 for an average 
background energy of 1.0 photon per detector element per slot. This Figure shows the symbol-error 
probability as a function of total average absorbed signal energy, K s , using estimates of the signal 
10 intensity obtained from the simulation with L = 1000, as well as three other curves where exact 
knowledge of the input distribution was assumed; the 16x16 subarray, the optimized subarray based 
on the actual energy probabilities, and the subarray obtained by maximizing SNRj. It should be 
noted that the simulation and SNRi subarray curves are indistinguishable in both cases, indicating 
that estimation of the signal energies over the array does not result in any appreciable performance 
15 degradation. It also is evident that subarray optimization based on the simple SNR i algorithm results 
M in only slight losses, but succeeds in greatly reducing the complexity of the estimator; for K b = 1.0, 
i|| the loss is less than 0.15 dB over the entire range considered in Figure 13. 

1^1 The elements of the apparatus and method and the general features of the components are 

Hi shown and described in relatively simplified and generally symbolic manner. Appropriate structural 
20 yj details and parameters for actual operation are available and known to those skilled in the art with 
j\ respect to the conventional aspects of the process. 

N Although specific embodiments are disclosed herein, it is expected that persons skilled in the 

|l| art can and will design alternative adaptive optical receivers that are within the scope of the 
Yl following claims either literally or under the Doctrine of Equivalents. 
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